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Abstract 

Superscaling analyses of inclusive electron scattering from nuclei are extended from the quasielas- 
tic processes to the delta excitation region. The calculations of (e, e') cross sections for the target 
nucleus ^■^C at various incident electron energies are performed using scaling functions f{ip') ob- 
tained in approaches going beyond the mean-field approximation, such as the coherent density 
fluctuation model (CDFM) and the one based on the light-front dynamics (LFD) method. The 
results are compared with those obtained using the relativistic Fermi gas (RFG) model and the 
extended RFG model (ERFG) . Our method utilizes in an equivalent way both basic nuclear quan- 
tities, density and momentum distributions, showing their role for the scaling and superscaling 
phenomena. The approach is extended to consider scaling function for medium and heavy nuclei 
with Z ^ N ioT which the proton and neutron densities are not similar. The asymmetry of the 
CDFM quasielastic scaling function is introduced, simulating in a phenomenological way the effects 
which violate the symmetry for if;' > including the role of the final-state interaction (FSI) . The 
superscaling properties of the electron scattering are used to predict charge-changing neutrino- 
nucleus cross sections at energies from 1 to 2 GeV. A comparison with the results of the ERFG 
model is made. The analyses make it possible to gain information about the nucleon correlation 
effects on both local density and nucleon momentum distributions. 

PACS numbers: 25.30.-c, 21.60.-n, 25.30.Pt, 21.10.Ft 
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I. INTRODUCTION 



Over the past four decades electron scattering has provided a wealth of information on 
nuclear structure and dynamics. Form factors and charge distributions have been extracted 
from elastic scattering data, while inelastic measurements have allowed for a systematic 
study of the dynamic response over a broad range of momentum and energy transfer. The 
nuclear i^scalin^ analysis of inclusive electron scattering from a large variety of nuclei (e.g. 

showed the existence of high-momentum components in the 
nucleon momentum distributions n{k) at momenta k > 2 fm~^ due to t he p resence of 

) that this 
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nucleon-nucleon (NN) correlations. It was shown (see, e.g. jlll 
specific feature of n{k), which is similar for all nuclei, is a physical reason for the scaling 
and superscaling phenomena in nuclei 



The concepts of scaling 



have been explored in 
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and superscaling 
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17| for extensive analyses of the (e, e') world data (see also |l8l|). 
Scaling of the first kind (no dependence on the momentum transfer) is reasonably good as 
expected, at excitation energies below the quasielastic (QE) peak, whereas scaling of second 
kind (no dependence on the mass number) is excellent in the same region. When both 
types of scaling behavior occur one says that superscaling takes place. At energies above 
the QE peak both scaling of the first and, to a lesser extent, of the second kind are shown 
to be violated because of important contributions introduced by effects beyond the impulse 
approximation, namely, inelastic scattering together with correlation contributions 

and meson exchange currents ^ I22 1 . 

The superscaling analyses of inclusive electron scattering from nuclei for relatively high 
energies (several hundred MeV to a few GeV) have recently been extended to include not 

23. A good 



only quasielastic processes, but also the region where A-excitation dominates 
representation of the electromagnetic response in both quasielastic and A regions has been 
obtained using the scaling ideas, importantly, with an asymmetric QE scaling function 
f'^^{ip') {4'' is the scaling variable in the QE region) and a scaling function f^{ip'/\) in the 
region up to inelasticities where the A contribution reaches its maximum. Both functions 
were deduced from phenomenological fits to electron scattering data. Particularly, for the 
scaling function in the quasielastic region it has been shown in Ref. |3] that, in contrast to 
the relativistic Fermi gas model scaling function, which is symmetric, limited strictly to the 
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region —l<'ip'<+l, and with a maximum value 3/4, the empirically determined f'^^ijp') 
has a somewhat asymmetric shape with a tail that extends towards positive values of ip' 
and its maximum is only about 0.6. Of course, the specific features of the scaling function 
should be accounted for by reliable microscopic calculations that take FSI into account. In 
particular, the asymmetric shape of f^^ tested in Refs. 2J, |25| by using a relativistic mean 
field (RMF) for the final states shows a very good agreement with the behavior presented 
by the experimental scaling function. 

The superscaling analyses and the present knowledge of inclusive electron scattering al- 
lowed one to start studies of neutrino scattering off nuclei on the same basis. The reactions 
of incident neutrino beams interacting with a complex nucleus have offered unique opportu- 
nities for exploring fundamental questions in different domains in physics. Recently, positive 
signals of neutrino oscillations confirmed the hypothesis of non-zero neutrino masses and trig- 
gered much interest on this issue 2^. To better analyze the next generation of high-precision 
neutrino oscillation experiments and to reduce their systematic uncertainty both neutral- 



[e.g. 
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I) and charged- current (e.g. 
neutrino-nucleus scattering have stimulated detailed investigations. 

The neutrino-nucleus interactions have been studied within several approaches investi- 
gating a variety of effects. Using the superscaling analysis of few-GeV inclusive electron- 
scattering data, a method was proposed in Ref. to predict the inclusive uA and VA 
cross sections for the case of ^^C in the nuclear resonance region, thereby effectively in- 
cluding delta isobar degrees of freedom. It was shown in Refs. that the important 
final-state interaction effects arising from the use of relativistic optical potentials within a 
relativistic Green's function approach lower the cross section by at least a 14% factor at 



incoming neutrino energies of 1 GeV. A similar result has been obtained in Refs. |3 
where the use of Random Phase Approximation (RPA) to predict the neutrino-nucleus cross 
section was discussed. Apart from relativistic dynamics and FSI, other effects may influ- 
ence the neutrino-nucleus reactions. The role of Pauli blocking and FSI in charged-current 

neut.„o.duced.ea.t,o,..a.a.y.ed,nRefs.0yr 

In this article we follow our method presented in Refs 
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to calculate the sca^ 



function in finite nuclei firstly within the coherent density fluctuation model (e.g., Refs. 



41 



mg 



40 



42L 43^). This approach, which is a natural extension of the RFG model, has shown how 



both basic quantities, density and momentum distributions, are responsible for the scaling 
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and superscaling phenomena in various nuclei. Although the scaling function obtained in 
[isj l is symmetrical around ip' = 0, the results agree with the available experimental data 
at different transferred momenta and energies below the quasielastic peak position, showing 
superscaling for ip' < including ip' < —1, whereas in the RFG model /{ip') = for 
■j/;' < — 1. It was shown in 3| that the QE scaling function can be obtained within the 
CDFM in two equivalent ways, on the basis of the local density distribution, as well as 
of the nucleon momentum distribution. As pointed out in [iJl, the nucleon momentum 

n 

distributions n(k) for various nuclei obtained in HSi within a parameter-free theoretical 

n n 

approach based on the light-front dynamics method (e.g., |4a] and references therein) 
can also be used to describe both y- and ■^/''-scaling data. So, in our present work we 
explore both methods, CDFM and LFD, to investigate further the scaling functions and 
their applications to analyses of electron- and neutrino scattering off nuclei. 

Our work is motivated by the fact that different models of the nuclear dynamics (in- 
cluding those with RMF dynamics and with RPA-type correlations accounted for) describe 
with different success the basic size and shape of the cross sections in studies of high-energy 
inclusive lepton scattering used so far. For this reason we extend further our consideration 
and calculate within the CDFM and LFD the scaling functions in the kinematical regions 
of the QE and A peak on the basis of momentum and density distributions of finite nu- 
clear systems in which nucleon correlations are included. This can be done either by using 
available empirical data for these quantities or theoretical calculations in which correlations 
are included to some extent. Then, the obtained scaling functions are applied to calculate 
electron-nucleus cross sections in QE and A regions in the energy range from 500 MeV to 
2 GeV for the target nucleus ^^C and to predict charge-changing neutrino and antineutrino 
reaction cross sections from the scaling region to the QE peak at energies of few GeV. We 
also make comparisons of the results obtained using our methods with those obtained using 
the RFG model and other theoretical schemes. 

The paper is organized in the following way. In Sec. II we present the formalism needed 
in studies of scaling functions in the quasielastic region and validate the superscaling within 
the CDFM and LFD for a variety of nuclei with Z = N and Z N. Then, we consider 
the nucleon momentum distributions and their applications in both approaches showing the 
sensitivity of the calculated scaling functions to the peculiarities of n{k) in different regions 
of momenta. Section III contains the CDFM and LFD methods to build up the scaling 
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function in the A region. The formahsm involved in obtaining the electron-nucleus cross 
sections in QE and A kinematical regions and the results of the numerical calculations are 
presented in Sec. IV A. In Sec. IV B we give our theoretical predictions for cross sections of 
quasielastic charge- changing neutrino reactions. Finally, in Sec. V we summarize the results 
of our work. 



II. SCALING FUNCTION IN THE QUASIELASTIC REGION 



A. QE scaling function in the CDFM 



As already mentioned in the Introduction, the superscaling behavior was firstly considered 



within the framework of the RFG model 



n, tne superscaling DeUi 

Q, H Q Q H Q 



where a properly defined 

function of the variable was introduced. As pointed out in Q], however, the actual 
nuclear dynamical content of the superscaling is more complex than that provided by the 
RFG model. It was observed that the experimental data have a superscaling behavior in the 
low-Co" side {uj being the transfer energy) of the quasielastic peak for large negative values 
of ip' (up to ip' ~ ~2), while the predictions of the RFG model are f{ip') = for ip' < — 1. 
This imposes the consideration of the superscaling in realistic finite systems. One of the 



approaches to do this was developed [13|, |l4| in the CDFM 



41 



42 



43| which is related to 
47| . It was shown in ^ , 14 , 3| 



the (5-function limit of the generator coordinate method 
that the superscaling in nuclei can be explained quantitatively on the basis of the similar 
behavior of the high-momentum components of the nucleon momentum distribution in light, 
medium and heavy nuclei. As already mentioned, the latter is related to the effects of the 



NN correlations in nuclei (see, e.g. 
The scaling i 



40 



4l|). 



The scaling function in the CDFM was obtained starting from that in the RFG model 
in two equivalent ways, on the basis of the local density distribution p(r) 
and of the nucleon momentum distribution n{k). This allows one to study simultaneously 
the role of the NN correlations included in p(r) and n{k) in the case of the superscaling 
phenomenon. To explore these properties the scaling function f{ip') has been derived in two 
ways in CDFM in [l^. Firstly, by means of the density distribution p{r), it leads to 

'■"/(fcFl^'l) 



/0^(^') 



dR\F{R)\'f^^^{^'{R)), 



(1) 
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with a weight function of the form 

\F{R)\' 

where 



1 dp{r) 



PoiR) 



Po{R) dr 
3A 



AttR^' 

fRFci'^'iR)) with ip'{R) = kpRip'/a is the scahng function related to the RFG model 

2" 



(2) 
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m^r being the nucleon mass and a = (977^/8)^/3 ^ 1.52^1/2. Secondly, by means of the 
momentum distribution n{k), the scaling function is expressed by 



/«-(^')= / dkF\G{kF)\'fSFGii^'{kp)), 



where ip'^kp) = kptp'/kp and the weight function is 

1 dn{k) 



(5) 



no{ki 



dk 



k=ki 
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(6) 



(7) 



In Eq. ® the RFG scaling function /?^g(V''(^f)) can be obtained from fSrci'^'iR)) [Eq- O] 
by changing there a/Rhj kp. In Eqs. (Q), and (0) the Fermi momentum kp is not a 
free parameter for different nuclei as it is in the RFG model, but kp is calculated within the 
CDFM for each nucleus using the corresponding expressions: 

1 o 47r(97r)i/3 



kp 

when the condition 
is fulfilled and 
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dRkp{R)\F{R)\'^ = a / dR^\F{R)\' 
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hm [piR)R'] = 
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dRp{R)R (8) 



(9) 
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dkpn{kp)kp 



(10) 



when the condition 
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kp^oo 



n{kF)kF^ 







(11) 



is fulfilled. 

As shown in the integration in Eqs. ((H) and using Eqs. ^ and (jHI), leads to the 
explicit relationships of the scaling functions with the density and momentum distributions: 



47r 
1 



a/(kpW\) 



dRp{R) 



R^dfZMM. 
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and 
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fQ^^^') = / dkMkF) 

^ Jkp\iP'\ 



kj^dj^g^igik^ 
3 dkp 



the latter at 



lim 

fci?— >oo - 



n{kF)kp^ 



:i2) 



(13) 



(14) 



One can see the symmetry in both Eqs. ()12|) and (fT^ written in r- and A;-space. We also 
note that in the consideration up to here the CDFM scaling function f^^{ip') is symmetric 
under the change of ip' by -%p' . 

In Refs. 3) Ql we used the charge density distributions to determine the weight function 
|F(i?)|2 and f^^{il)') in Eqs. (H)), and © for the cases of ^He, ^^C, ^^Al, sspe and ^^^Au. 
The results for the scaling function f^^i^ip') agree well with the available data from the 
inclusive quasielastic electron scattering for "^He, ^^C, ^^Al, ^^Fe and only approximately 
for ^^"^Au for various values of the transfer momentum q = 500, 1000, 1650 MeV/c |l3| and 

, showing superscaling for negative values of ip' including also those 
smaller than -1, whereas in the RFG model f{ip') = for t/^' < — 1. One can see this in 
Fig. ID for ^He, and ^^Al at g = 1000 MeV/c. At the same time, however, in 0Q we 
encountered some difficulties to describe the superscaling in ^^^Au which was the heaviest 
nucleus considered. We related this in to the particular A-dependence of n{k) in the 

model that does not lead to realistic high-momentum components of n{k) in the heaviest 



nuclei. We followed in Refs 



.0 



14| an artificial way to "improve" the high- momentum tail 



of n{k) in ^^''Au by taking the value of the diffuseness parameter b in the Fermi-type charge 
density distribution of this nucleus to be 6=1 fm instead of the value 6=0.449 fm (as obtained 

n 

from electron elastic scattering experiments, see e.g. j48[). In this way the high- momentum 
tail of n{k) for ^^^Au in CDFM becomes similar to those of ^He, ^^C, ^^Al, and ^^Fe and this 
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leads to a good agreement of the scaling function /'^^{ip') with the data also for ^^^Au. Still 
in we pointed out, however, that all the nucleons (not just the protons) may contribute 
to f^^iip') for the transverse electron scattering and this could be simulated by increasing of 
the diffuseness of the matter density with respect to that of the charge density for a nucleus 
like ^^^Au that has much larger number of neutron than of protons. 
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FIG. 1: The quasielastic scaling function f^^{ilj') at q = 1000 MeV/c for ^He, ^^c, 27^1, 82j^r 
and ^^■^Au calculated in CDFM. Dotted line: RFG model result. The curves for ^He, and ^^Al 
nuclei almost coincide. Grey area: experimental data 
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In jl^ we assumed that the reason why the CDFM does not work properly in the case of 



^^''Au is that we had used in 



14| only the phenomenological charge density, while this 
nucleus has many more neutrons than protons (A^=118 and Z=79) and therefore, proton 
and neutron densities may differ considerably. In the case when Z ^ N and the proton and 
neutron densities are not similar, the total scaling function will be expressed by the sum 
of the proton fp^{ip') and neutron f^^iip') scaling functions which are determined by the 
proton and neutron densities Pp{r) and Pn{r), respectively: 



«p(n)/(fc?.'"'l^'l) 



41 (^') 



(15) 



In Eq. ()15j) the proton and neutron weight functions are obtained from the corresponding 
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proton and neutron densities 

,2 ArcR^ dpp^n){r) 



\Fp{n){,R)\ 



3Z(iV) dr 



(16) 



r=R 



oo 

J p,(„)(r)cir = Z(iV), (18) 



and the Fermi momentum for the protons and neutrons is given by 

oo 

fc^^"^ = «p(n) j (19) 


The RFG proton and neutron scahng functions f^p]j{ip' {R)) have the form of Eq. where 
a and kp stand for from Eq. (fTTj) and fc^*-""-* from Eq. (fTUj) . respectively. The functions 
are normahzed as follows: 

oo 

j \Fp^^){R)\'dR=l, (20) 



oo 

eQE 



/i)(^')^^'=l- (21) 

Then the total scaling function can be expressed by means of both proton and neutron 
scaling functions: 

/Q^(V^') = \[Zf?^{i;') + iV/f ^(^')] (22) 

and is normalized to unity. 

The same consideration can be performed equivalently on the basis of the nucleon mo- 
mentum distributions for protons nP{k) and neutrons n"'{k) presenting /^^{ip') by the sum 
of proton and neutron scaling functions ()22|1 calculated similarly to Eqs. p5|l - (f^ (and to 
Eqs. ©, ©, (HDD and (HH)): 

oo 

/ph(^')= / rffcF|Gp(.)(:fci.)|V^?i(^'(^F)), (23) 

where _ 

. ^ (24) 



3Z{N) dk 
10 



k=kp 



with f^Fci'ip'ikF)) containing ap(„) from Eq. (fTTj) and fc^*-"^ calculated as 

oo 

}4^) = J dkFkF\Gp(^n){kF)?. (25) 


The scaling functions for several examples, such as the medium stable nuclei ^^Ni and ^^Kr 
and the heavy nuclei ^^^Sn and ^^^Au are calculated following Eqs. ()15|) - (P^ using the 
corresponding proton and neutron densities obtained in deformed self-consistent mean-field 
(HF+BCS) calculations with density-dependent Skyrme effective interaction (SG2) and a 
large harmonic-oscillator basis with 11 major shells 0, Is^. In Fig. ^ we give the results 
for the ^^Kr and ^^^Au nuclei in which Z ^ N and compare them with the results for 
^He, ^^C {Z = N) and ^'^Al {Z ~ TV), as well as with the experimental data (presented 
by a grey area and taken from Q) obtained for ^He, ^^C, ^^Al, ^<^Fe, and '^''An. The 
scaling functions are in a reasonable agreement with the data, which was not the case for 
^^''Au calculated in by using only the Fermi-type char ge d ensity with phenomenological 
parameter values b = 0.449 fm and R = 6.419 fm from j48|. At the same time we note 
also the improvement in comparison with the RFG model result in which f'^^{ip') = for 
ip' < —1. Thus, it can be concluded that the scaling function f^^{'4)') for nuclei with Z ^ N 
for which the proton and neutron densities are not similar has to be expressed by the sum 
of the proton and neutron scaling functions. The latter can be calculated by means of 
theoretically and/or experimentally obtained proton and neutron local density distributions 
or momentum distributions. 

As known (e.g. jl^ . Ssl), the total inclusive electron scattering response is assumed to be 
composed of several contributions: i) the entire longitudinal contribution which superscales 
and is represented by the QE scaling function f^^{ip')'-, ii) a part of the transverse response, 
which arises from the quasielastic knockout of nucleons and is also driven by the scaling 
function f^^{ip'), and iii) the additional contribution of the transverse response from MEG 
effects and from inelastic single-nucleon processes including the excitation of the A isobar. 
The effects of point iii) break the scaling. In an universal scaling function f^^{ip') has 
been determined by reliable separations of the empirical data into their longitudinal and 
transverse contributions for A > A. Such separations are available only for a few nuclei 
[sij ]. All of these response functions have been used to extract the "universal" QE response 
function /^^{ip') (see Fig. 1 of [2^) which is parametrized by a simple function. This 



11 



function has a somewhat asymmetric shape. Its left tail {ip' < 0) passes through the grey 
area of Fig. [TJ The right tail {ip' > 0) extends larger towards positive values of ip'. In 
contrast, the RFG scaling function is symmetric. The CDFM scaling function discussed so 
far, which is based on the RFG one, is also symmetric. As mentioned, the maximum value 
of f^^iip') in RFC (and in CDFM) is 3/4, while the empirical scaling function extracted in 
Ref. j23j reaches about 0.6. 
As mentioned in t 



and relativized 



, if FSI are neglected, the RMF theory 
shell-model studies |54] provide rather modest differences from the RFG predictions. An- 
other possible reason for the differences between the RFG (or mean-field results) and the 
empirically determined scaling function arises from high-momentum components in realistic 
wave functions which may be large enough. In Q the scaling function was taken from the 
experiment. In the present work we also limit our approach to phenomenology when con- 
sidering the asymmetric shape and the maximum value of the quasielastic scaling function. 
In order to simulate the role of all the effects which lead to asymmetry, we impose the latter 
on the RFG scaling function (and, correspondingly, on the CDFM one) by introducing a 
parameter which gives the correct maximum value of the scaling function (ci in our expres- 
sions given below) and also an asymmetry in f^^{ip') for tp' > We consider the main 
parts of the RFG scaling function for ip' <Q and ip' > in the following forms, keeping the 
parabolic dependence on ip' as required in Ref. |lO|: 



(26) 



The total RFG scaling function is normalized to unity: 



\p-&\ 


e 


1 - 


0"1 











ip' > 0. 



cQE 



If the normalization of the scaling function for negative values of ip' is equal to 



2 

di^'fRFcM') = qCl, 

-: . O 



(27) 



(28) 



(29) 



then, in order to keep the total normalization [Eq. (j28|) ]. the normalization for positive ip' 
has to be: 

l-a = y^ rf^7?FG,2(^') = 3C1C2. (30) 
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From Eqs. ()29|) and ()3U|) we get the relationship between C2 and Ci: 

2ci 



(31) 



In the RFG ci = 3/4 and, correspondingly, C2 = 1. In the CDFM the QE scaling function 
will be: 

f'^''{^') = f?''{^') + f?\^'), (32) 

where 



i'a/kp\il)'\ 

/f^(^')= / dR\F{R)\'ci 



a 



C2a/kp\il)'\ 



dR\F{R)\^ci 



1 - 



kFRW\ 
C2a 



(33) 



(34) 



In this approach, parametrizing the RFG scaling function by the coefficient ci we account 
for the experimental fact that Ci ^ 3/4 and take this value in accordance with the empirical 
data. Then from the normalization [Eqs. (j28p - (j3(J|) ] we determine the corresponding value of 
C2 using Eq. dnH). As in |l3, 14|, the CDFM scaling function is obtained [Eqs. ((221) -(El] by 
averaging the RFG scaling function. As an example, we give in Fig. |21the CDFM QE scaling 
function for different values of Ci (0.75, 0.72, 0.60 and 0.50) in comparison with the empirical 
data and the phenomenological fit. We also include for reference the scaling function ob- 



tained from calculations for (e, e') reaction based on the re^ 



n 



ativistic impulse approximation 
25| for details). In this way we 



(RIA) with FSI described using the RMF potential (see 
simulate in a phenomenological way the role of the effects which violate the symmetry for 
positive values of ip' of the QE scaling function, which in the RMF approximation are seen 
to be due to the FSI. 



B. QE scaling function in the LFD method 

In this Subsection we will obtain the QE scaling function on the basis of calculations of 
nucleon momentum distribution (using Eqs. (5)-(7) or Eq. (13)) obtained within a modi- 
fication of the approach from 441. The latter is based on the momentum distribution in 

n n 

the deuteron from the light-front dynamics (LFD) method (e.g., |4^ |4^ and references 
therein). Using the natural-orbital representation of the one- 
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FIG. 2: The quasielastic scaling function f'^^i'tp') for calculated in CDFM m comparison 



with the experimental data (black squares) 



23]. The CDFM results for different values of ci are 



presented by solid lines. Also shown for comparison is the phenomenological curve which fits the 
data (dash-two dots), as well as the curve (dash-dot line) corresponding to the result for (e, e') 
obtained within the relativistic impulse approximation and FSI using the relativistic mean field 

nn 

(see Refs. |2j,|25|). 

was written as a sum of contributions from hole-states [n'^(/i;)] and particle-states [nP(/c)] 
(see also Q) 

riAik) = Na [n''{k) + nP(fc)] . (35) 

In dSH) 

F.L. 

n^A;) = C(A;) 5^2(2j + l)A„,,|/?„i,(A;)p, (36) 

nlj 

where F.L. denotes the Fermi level, and 



C{k) 



rriN 



2 ' 



(37) 



rriN being the nucleon mass. To a good approximation for the hole states, the natural 
occupation numbers Xnij are close to unity in j4J| and the natural orbitals Rnij{k) are 
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replaced by single-particle wave functions from the self-consistent mean-field approximation. 
In [4^] Woods-Saxon single-particle wave functions were used for protons and neutrons. A^^ 
is the normalization factor. Concerning the particle-state [nP(/i;)] contribution in (j35p . we 
used in [44] and the well-known facts that: (i) the high-momentum components of 
n{k) caused by short-range and tensor correlations are almost completely determined by the 
contributions of the particle-state natural orbitals (e.g. (ssl), and (ii) the high- momentum 
tails of nA{k)/A are approximately equal for all nuclei and are a rescaled version of the 
nucleon momentum distribution in the deuteron nd{k) 

UAik) ~ aAUdik), (38) 



where a a is a constant. These facts made it possible to assume in 



44 and H that n^ik) 



is related to the high-momentum components n^{k) of the deuteron, that is, 

n^k) = ^n,{k). (39) 



In dSni) n^lk) is expressed by an angle-averaged function [44] as 



n,ik)=C{k){l-z^)fi{k). (40) 

In (j4Up z = cos(/c, n), n being a unit vector along the three vector (c^) component of the four- 
vector u which determines the position of the light-front surface 45|,|46|. The function f^lk) 
is one of the six scalar functions fi-eik'^ ^n^k) which are the components of the deuteron total 
wave function \E'(fc,n). It was shown [45] that /s largely exceeds other /-components for 
k > 2.0-2.5 fm~^ and is the main contribution to the high-momentum component of nd{k), 
incorporating the main part of the short-range properties of the nucleon-nucleon interaction. 

It was shown in Fig. 2 of that the calculated LFD n{k)^s are in good agreement with 
the "|/-scahng data" for "^He, ^^C and ^^Fe from P] and also with the yew analysis 0,0] up 
to A; < 2.8 fm~^. For larger k the momentum distributions from LFD exceeds that obtained 
from yew analysis. We should note also that the calculated scaling function f^^i^^j') using 
the approximate relationship (see Eq. (75) and Fig. 4 of [3]) 

J\y\ 2c 2mM 

for ^^Fe at g = 1000 MeV/c is in agreement with the data for —0.5 < V^' < 0, while in the 
region —1.1 < ip' < —0.5 it shows a dip in the interval —0.9 < ip' < —0.6. This difference 
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is due to the particular form of n{k) from LFD shown in Fig. 2 of [14] (a dip around 
k ^ 1.7 fm~^ and a very high-momentum tail at > 2.8 fm~^). This result showed that the 
assumption (j39|) for the particle-state contribution is a rather rough one. In this paper we 
consider a modification of the approach in which we include partially in the particle-state 
part n^lk) not only n^^lk) but also n2{k) which is related to the angle-averaged function 

f2{ky. 



n2{k)=C{k)fl{k). (42) 
Then the particle-state part can be written in the form 

n''{k) = P[n2{k)+n,{k)], (43) 

where /3 is a parameter. Then the LFD nucleon momentum distribution for the nucleus with 
A nucleons will be: 

ULFDik) = NA[n''{k)+(3{n2{k)+n,{k))], (44) 
with n^{k) from Eq. (j36|) and 



NA = Un dqq 







■F.L. n ^ -1 

^2(2j + l)Xni,C{q)\RniM\'+P{Ml)+MQ)) 

■ nlj 



(45) 



In Fig. 01 we present the nucleon momentum distribution for ^^C calculated within the LFD 
method using Eqs. §^-§7i), with the parameter value /? = 0.80. It is 



compared with the band of CDFM momentum distri 



(grey area), with ncw(^) from the yew analysis 



m 



butions for ^He, ^^C, 27A1, sepe, 197^^ 



and with the y-scaling data l3\ for 



4 



He, ^^C, and ^^Fe. It can be seen that up to A; ^ 2.8 fm ^ ^lfd curve is close to the results 



of . For the region 1 < k < 2.5 fm ^ it is between them and 



br k ^ 2.8 fm ^ it is 



close to ncw{k), in contrast to our previous results in Fig. 2 of Ref. Q (see also H) 

which 

were based on Eq. ()39|) and which are also shown for comparison in Fig. El This behavior 
of nLFD(^) reflects in the result of the calculation of the QE scaling function using Eq. pij) 
given in Fig. HI It can be seen that both momentum distributions new and ni^FD{k) 
[Eq. ()44|) ] give a good agreement with the experimental data for the QE scaling function at 
least up to ip' ~ —1.2. This result is an improvement of that for LFD shown in Fig. 4 of jl^ . 
where only the contribution was used in the calculation of nP{k) and nLFD(/^)- 
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FIG. 3: The nucleoli momentum distribution n{k). Grey area: CDFM combined results for ^He, 
^■^C, ^^Al, ^^Fe and ^^^Au. Solid line: result of the present work for ^■^C using the modified LFD 
approach with /3 = 0.80. Dashed line: yew-scaling result [t, Dash-dotted line: result of 
LFD for ^^C from Dotted line: the mean-field result using Wood-Saxon single-particle wave 



functions for """Fe. Open squares, circles and triangles are y-scaling data [5|] for He, C and "'"Fe, 
respectively. The normalization is: / n(A;)(i^k = 1. 

III. SCALING FUNCTION IN THE QUASIELASTIC DELTA REGION 



In this Section we will extend our analysis within both CDFM and LFD to the A-peak 
region, which is not too far above the QE peak region and is the main contribution to 
the inelastic scattering. Dividing the cross section by the appropriate single-nucleon cross 
section, now for N —>■ A transition, and displaying the results versus a new scaling variable 
(iJj 'a) ( in which the kinematics of resonance electro-production is accounted for) it is obtained 
in j23| that the results scale quite well. This is considered as an indication that the procedure 
has identified the dominant contributions not only in the QE region, but also in the A-region. 
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FIG. 4: The quasielastic scaling function f'^^W) calculated using Eq. (|IT|) at q = 1000 MeV/c 
with ncw{k) from the yciy-scaling analysis 0,0] for ^^Fe (solid line) and niFDik) from modified 
LFD approach [Eq. (jlH)] for ^^C (dashed line). 

The shifted dimensionless scahng variable in the A-region ip'^ is introduced (see, e.g. [2^) 
by the expression: 



where 



.1/2 



+ 1, A' > 
-1, A' < A' 



rriN 



A' = A- 



Eshift 



t' = k^- A'2 



(46) 



(47) 



(48) 



A 



2m 



N 



2m 



2 \ 2 
T = K — X 



N 



hift 



2m N 



-^A 



^l\ + 4/€2 - 1 



(49) 
(50) 



/?A 



PA = 1 + — , Pa = 1 + — , 
r r' 



(51) 
(52) 
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/3a = ^ (/ii - 1) . (53) 
The relativistic Fermi gas superscaling function in the A domain is given by 

/HFG(^k) = ^(l-V^k')0(l-V^k')- (54) 



Following the CDFM application to the scaling phenomenon, the A-scaling function in 
the model will be: 



/^(^k)= / dR\F^{R)\'fi^^{f^{R)). 



In Eq. dSni): 



1 + 



m 



N 



f^yPA + — - A Pa - 1 



where 



- 1 



m 



N 



a 



'1 + 



kUR) 



m 



N 



and kpiR) = —■ 
R 



(55) 



(56) 



(57) 



In the CDFM kp can be calculated using the density distribution (Eqs. ((HI), © or 
and (HSl)) or the momentum distribution (Eqs. (fTUj). (fTTj) or and (jHI)). The weight 
function |Fa(-R)P is related to the density distributions (Eqs. Q or |T6|) ). In the equivalent 
form of the CDFM, the scaling function can be written in the form: 



dkF\G^{kF)\^f^pa{MkF)) 



where G/^lkp) is determined by means of the momentum distribution (Eqs. © or 



(58) 



and 



i;'l{kp) = t{kF).iJ. 



12 



(59) 



with 



t{kF) = 



m 



N 



(60) 



m 



N 
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Here we would like to note that though the functional forms of /^("^Aa) [Eq- ()55|) ] and the 
weight function |Fa(-R)P (Eqs. (0) or (fTBj) ) are like before, i.e. as in the case of the QE region, 
the parameters of the densities (e.g. the half-radius and the diffuseness when Fermi- 
type forms have been used) may be different from those {R and b) in the QE case. Along 
this line, we calculated firstly the scaling function f^^ip'^) by means of Eqs. (fH ^ -(|37 |) using 
the Fermi-type density for ^^C. We found the values of Ra and 6a fitting the scaling data 
at the A peak extracted from the high-quality world data for inclusive electron scattering 
(given in Fig. 2 of Our results are presented in Fig. El As mentioned already in the 

QE case, the empirical data require to use a value of the coefficient in the right-hand side 
of Eq. for the REG scaling functions f^pci'^'A) different from 3/4. In our calculations 
in the A-region we use the value 0.54. We found that reasonable agreement with the data 
can be achieved using the parameter values Ra = 1.565 fm and 6a = 0.420 fm (the Fermi- 
momentum value is taken to he kp = 1.20 fm~^ and this choice leads to normalization 
to unity of fRpci'^'A))- The value of -Ra is smaller than that used in the description of 



the QE superscaling function for ^^C [13, {R = 2.470 fm) while the value of 6a is the 



same as 6 in the QE case. Secondly, we calculated /^(V'a) using Eqs. (jKH|) . (fHT?|l and ^K^ . 
In Eq. ()58|) the weight function |G'a(^f)P was determined by means of Eq. © and the 
nucleon momentum distribution nLFD (Eqs. (jl^ and (US))) calculated with the parameter 
value P = 0.80 (shown in Fig. E)). We note that the use of ni^FD{k) with this value of /3 gives 
simultaneously a reasonable agreement both with the results for the momentum distribution 
from the y-scaling data shown in Fig. |21 as well as with the QE scaling function shown in 



IV. SCALING FUNCTIONS AND INCLUSIVE LEPTON SCATTERING 
A. Scaling functions and (e, e') reaction cross sections 

In the beginning of this Subsection we will give some basic relationships concerning 
inclusive electron scattering from nuclei. An electron with four-momentum k^^ = (e, k) is 
scattered through an angle 6 to four-momentum k'^ = (e', k'). The four-momentum transfer 
is then 




Fig.ia 



Q'' = ik-kr = iu;,q), 



(61) 
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FIG. 5: The f^i^p'^) scaling function for in the A -region. Dashed hne: CDFM result (with 
i?A = 1.565 fm, 6a = 0.420 fm, kp = 1.20 fm~^). Solid line: result of modified LFD approach 
(f3 = 0.80, kp = 1.20 fm~i). The coefficient ci^O.54 in both CDFM and LFD cases. Averaged 
experimental values of /^(V'a) taken from 



where u = e — e', q = |q| = |k — k'| and 



(62) 



In the one-photon-exchange approximation, the double-differential cross section in the 
laboratory system can be written in the form (e.g. jlo|): 



— ) Riiq,^) 



2 2 ' ^'^^^'^> 



where 



acos(e/2) 



(63) 



(64) 



.2esin2(^/2). 

is the Mott cross section and a is the fine structure constant. In Eq. ()63p Rp and Rt are 
the longitudinal and transverse response functions which contain all the information on the 
distribution of the nuclear electromagnetic charge and current densities, being projections 
(with respect to the momentum transfer direction) of the nuclear currents. They can be 
separated experimentally by plotting the cross section against tan^(6'/2) at fixed (g, cu) (the 
so-called "Rosenbluth plot"). These functions can be evaluated as components of the nuclear 
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tensor W^i,. In [l^ this tensor is computed in the framework of the RFG model and Rl{t) 
for the QE electron scattering are expressed by means of the RFG scaling function (Eq. (jH)) 
of Ref. [lo|). 

At leading-order in the parameter kp/niiq the QE responses have the form 



A) = Ao^[(l + T)W,{r) - W,{r)] x /^^^ (65) 



A) = Ao[2iyi(T)] X /g^clV''), (66) 

with 

Ao = 3-, (67) 

where M = Z 01 N and Wi^ W2 are the structure functions for elastic scattering which are 
linked to the Sachs form factors 

(1 + t)W2{t) - W^{t) = Glir), (68) 



2Wi{t) = 2tGI,{t). (69) 

In 00 

the electro-production of the A-resonance is considered computing the nuclear 
tensor also within the RFG model and analytical expressions for the response functions are 
obtained. The latter contain the RFG A-peak scaling function (I54j) and read 0|: 

Rd^A) = ^^^^[il + Tp')w2iT)-w^{T)+W2iT)DiK,\)] X fipaH^'^), (70) 

Rt{^, A) = /-^^^ [2w,{r) + W2{r)D{K, A)] x /^^^(V^k), (71) 
where Af = Z or N, 

D{k, A) ^ ^[(Ap + 1)2 + (Ap + 1)(1 + (72) 



The single-baryon structure functions can be expressed by means of the electric {Ge), mag- 
netic (Gm) and Coulomb (Gc) delta form factors 0|: 

wi(r) = ^(pA + l)'(2rp + 1 - pa)(G'1/ + 3G|), (73) 
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Mr) = + ^r^^ 

These form factors are parametrized as follows 



r 



(74) 



GMiQ') = 2.97 f{Q' 



(75) 



GEiQ') = -0.03/(g2), (76) 
G'c(Q') = -O.lSGAf (Q'), (77) 

where 



/(Q^) = G^^iQ')- ^ (78) 



3.5(GeV/c) 



with 



. - (1 + 4.97r)2 
being the Galster parametrization [58] of the electric form factor. 

In the CDFM the longitudinal and transverse response functions can be obtained by 
averaging the RFG response functions in the QE region [Eqs. and and A-region 
[Eqs. flTOj) and (f7T|) ] by means of the weight functions in r-space |F(i?) p and fc-space 1(7(^^7) p, 
similarly as in the case of the QE- and A-scaling functions (Eqs. (P), (jSJ, (fT3j) . (PHj) and (j^ . 
(j58|) . respectively). As a result, accounting for the different behavior of the RFG scaling 
functions and terms containing rjp{R) = kp{R)/mj^ as functions of R or kp = a/R in 
(jnSl), (inni), and (IHI), the CDFM response functions Rl(t) in QE- and A-regions have 
approximately the same forms as in the equations just mentioned, in which, however, the 
RFG scaling functions are changed by the CDFM scaling functions obtained in Sections 2 
and 3. 

In Figs. lUlfT^l we give results of calculations within the CDFM of inclusive electron scat- 
tering on ^^C at different incident energies and angles. The QE-contribution is calculated 
using the Fermi-type density distribution of ^^C with the same values of the parameters as 
in Q|: R = 2.47 fm and b = 0.42 fm (which lead to a charge rms radius equal to the 
experimental one) and Fermi momentum kp = 1.156 fm^^. The delta-contribution is calcu- 
lated using the necessary changes of the parameter values of the Fermi-type density (used 
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in Fig. El): -Ra = 1-565 fm, &a = 0.42 fm and kp = 1-20 fm~^. The coefficient Ci used in the 
A-region scahng function is fixed to be equal to 0.54 so that the maximum of the scahng 
function to be in agreement with the data. The scahng function /^(V'^) is symmetric, its 
maximum is chosen to be 0.54 (but not 0.75) and it is normahzed to unity by means of the 
fixed value of kp = 1.20 fm~^. The inclusive electron-^^C scattering cross sections shown in 
Figs. IBHT^ are the sum of the QE and A-contribution. The results of the CDFM calculations 
are presented for two values of the coefficient ci in the QE case (noted further by c^^), 
namely for c^^ ^ 0.72 and c^^ = 0.63. This is related to two types of experimental data. 
In the first one the transferred momentum in the position of the maximum of the QE peak 
extracted from data (uS^f) is > 450 MeV/c ~ 2kp, roughly corresponding to the domain 



where scaling is fulfilled |l9|, |23|. Such cases are presented in Figs. lUHT^ In these cases we 
found by fitting to the maximum of the QE peak the value of cf^ to be 0.72-0.73, i.e. it is 
not the same as in the RFG model case (case of symmetry of the RFG and of the CDFM 
scahng functions with cf'^ = 0.75), but is slightly lower. This leads to a weak asymmetry 
of the CDFM scaling function for cases in which > 450 MeV/c. In the second type of 
experimental data is not in the scaling region (g£.p < 450 MeV/c). Such cases are given 
in Figs. IT^lfTSl For them we found by fitting to the maximum of the QE peak extracted 
from data the value of cf^ to be 0.63. For these cases the CDFM scaling function is def- 
initely asymmetric. So, the results in Figs. IU1IT31 are presented for both almost symmetric 
(cf^ ~ 0.72) and asymmetric {cf^ = 0.63) CDFM scaling functions. One can see that the 
results for the almost symmetric CDFM scaling function agree with the electron data in the 
region close to the QE peak in cases where > 450 MeV/c and overestimate the data for 
cases where approximately qf^^ < 450 MeV/c. The results with asymmetric CDFM scaling 
function agree with the data in cases where < 450 MeV/c and underestimate the data 
in cases where > 450 MeV/c. Here we would like to emphasize that, in our opinion, 
the usage of asymmetric CDFM scaling function is preferable, though the results in some 
cases can underestimate the empirical data, because other additional effects, apart from QE 
and A-resonance (e.g. meson exchange currents effects) could give important contributions 
to the cross section for some specific kinematics and minor for others. A similar situation 
occurs for the results obtained within the RMF approach j25|| particularly when the CC2 
current operator is selected. 

In Table E we list the energies, the angles, the values of cf ^ obtained by fitting the 
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magnitude of the QE peak, and the energy shifts in the QE and A-case, as well as the 
approximate values of the transfer momentum g^.^ in the position of the maximum of the 
QE peak (cu^p) for different cases. The values of the energy shifts e^^f^f ■* for the QE- and 
A-regions are generally between 20 and 30 MeV. In the Figures we also present the QE- 
contribution (as well as A-contribution) for the value of c^^ which fits approximately the 
magnitude of the QE peak. 

2.0 



> 

^ 1.5h 



CO 



1.0 - 



3 
■a 

§ 0.5 h 
■a 



0.0 



-i 1 1 r 




200 400 600 

CD [MeV] 



800 



FIG. 6: Inclusive electron scattering on at e = 1299 MeV and 9 = 37.5° = 792 MeV/c 

> 2kp). The results obtained using c^^ = 0.72 in the CDFM scaling function for the QE cross 
section and the total result are given by dashed and thick solid line, respectively. Dotted line: 
using CDFM A-scaling function; thin solid line: total CDFM result with = 0.63. Dash-dotted 



line: result of ERFG method 



IS 



23l |. The experimental data are taken from 



In Figs. IHl and ^2 we present also the calculations of the electron cross sections using QE- 
and A-scaling functions obtained by using the nucleon momentum distributions obtained 
in the LFD approach (Section 3) which give a reasonable agreement with the empirical 
electron scattering data. In Figs. Inland El we also give for comparison the results of the 
cross sections obtained within the ERFG method jiol . [2^. In this method the response 
functions and differential cross sections are calculated using the scaling function fitted to 
the data. 

It is interesting to note that for those kinematics where the overlap between the QE and 
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FIG. 7: Inclusive electron scattering on at e = 2020 MeV and 6 = 20.02° (g£f = 703 MeV/c 
> 2kp). The results obtained using c^^ = 0.73 in the CDFM scaling function for the QE cross 
section and the total result are given by dashed and thick solid line, respectively. Dotted line: using 
CDFM A-scaling function. Thin solid line: total CDFM result with = 0.63. The experimental 



data are taken from 



6C|. 



A peaks is bigger (Figs. IHl [Z|and|H)), the asymmetric CDFM model {cf^ = 0.63) gives rise 
to an excess of strength in the transition region. This makes a difference with the ERFG 
model (see Fig.lHl) which fits nicely the data in that region. This discrepancy between the two 
models, asymmetric CDFM and ERFG, can be explained by noting the different behavior 
presented by the two scaling functions in the region of ip' between 0.5 and 1.5, being the 
asymmetric CDFM one significantly larger. 

Note on the other hand, that in the cases where the overlap between the QE and A peaks 
is weaker (Figs. lU HT^ . the asymmetric CDFM model, compared to the almost symmetric 
CDFM one, reproduces better the data in the transition region although it underpredicts 
importantly the maximum of the QE peak. Concerning results in Figs. ITI^llTHl (it can be also 
applied to Figs. fTD and fT^ . one observes that both CDFM approaches do not reproduce 
the strength of data located in the region between the QE and delta peaks. This is not the 
case for the ERFG model (see Fig. IT^ which fits nicely the experiment for uj > 180 MeV. 
This result is connected with the much bigger strength of the scaling function provided by 
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FIG. 8: The same as in Fig.Elfor e = 1108 MeV and 9 = 37.5° (g£p = 675 MeV/c > 2kF). Dot- 



dashed Hne: using Q 
data are taken from 



3- 



and A-scaling functions obtained in the LFD approach. The experimental 



the ERFG model for larger values of the scaling variable, ip' > 2 (see Fig. 

From this whole analysis, one may conclude that the phenomenological procedure intro- 
duced in the CDFM model to get an asymmetric scaling function, gives rise to an excess of 
strength in the region 0.5 < ip' < 1.5, whereas the model lacks strength for larger ip'-values, 
ip' > 2. 



B. Scaling functions and charge-changing neutrino-nucleus reaction cross sections 



In this Subsection we will present applications of the CDFM and LFD scaling functions 
to calculations of charge- changing neutrino-nucleus reaction cross sections. We follow the 
description of the formalism given in 23[ . The charge-changing neutrino cross section in the 
target laboratory frame is given in the form 



dQdk' 



X' 



(80) 
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FIG. 9: The same as in Fig. [71 for e = 620 MeV and 6 = 60° = 552 MeV/c > 2kF). The 



experimental data are taken from 
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FIG. 10: The same as in Fig. El for e = 2020 MeV and 6 = 15.02° = 530 MeV/c > 2^^) for 



the CDFM results. The experimental data are taken from 



6C| 



where x = + for neutrino- induced reaction (e.g., vi + n ^ + p, where 
X = — for antineutrino- induced reactions (e.g., + p + n), 



(70 



27r2 



■[A;'cos^/2]^ 



e, /i, r) and 



(81) 
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FIG. 11: The same as in Fig.Elfor e = 500 MeV and 6 = 60° (g?^ = 450 MeV/c > Ikp). Here the 
dot-dashed line shows the result using QE- and A-scahng functions obtained in the LFD approach. 
The experimental data are taken from j62l |. 

G = 1.16639 X 10"^ GeV~^ being the Fermi constant, 6c being Cabibbo angle (cos^c = 
0.9741), 



tan^ 9/2 = 



\Ql 



(82) 



^^0 



e + e')'-g2 = 4ee'-|g|^ 



(83) 



The quantity JF^ which depends on the nuclear structure is written in j23| as a gen- 
eralized Rosenbluth decomposition having charge-charge, charge-longitudinal, longitudinal- 
longitudinal and two types of transverse responses. The nuclear response functions are 
expressed in terms of the nuclear tensor W'^'^ in both QE and A-regions using its rela- 
tionships with the RFG model scaling functions. Following 2311 . in the calculations of the 
neutrino-nucleus cross sections the Hohler parametrization 8.2 j65| of the form factors in the 
vector sector was used, while in the axial-vector sector the form factors given in [2^ were 
used. 

In our work, instead of the RFG scaling functions in QE- and A-regions, we use those 
obtained in the CDFM and LFD approach (Sections 2 and 3). In Fig. we give the 
results of calculations for cross sections of QE neutrino (i/^,/x~) scattering (Figs. IT^ . c, 
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FIG. 12: The same as in Fig.Elfor e = 730 MeV and 9 ^37.1° {q9xy = 442 MeV/c < 2kF) for the 
CDFM results. The experimental data are taken from 
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FIG. 13: The same as in Fig. El for e = 1650 MeV and 9 = 13.5° (g£p = 390 MeV/c < 2kF) for 



the CDFM results. The experimental data are taken from 



64|. 



d, e, f) on ^^C and also antineutrino (77^, scattering fFig. 116b) for energies of neutrino 
= 1,1.5 and 2 GeV and of antineutrino = 1 GeV. The presented cross sections are 
functions of muon kinetic energy. The energy shift is equal to 20 MeV. The calculations of 
the neutrino-nucleus cross sections in the A-region will be a subject of a future work. 
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FIG. 14: The same as in Fig-Hfor e = 1500 MeV and 9 = 13.5° (g^p = 352 MeV/c < Ikp). The 
experimental data are taken from 64 1. 
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FIG. 15: The same as in Fig.Elfor e = 537 MeV and 9 = 37.1° (g£p = 326 MeV/c < 2kF) for the 
CDFM results. The experimental data are taken from {5^ . 

We give the results of our calculations using the CDFM scaling function which is almost 
symmetric (with ci = 0.72), as well as the asymmetric CDFM scaling function (with ci = 
0.63). These values of ci correspond to the cases of inclusive electron scattering considered. 
As can be seen the results obtained by using the almost symmetric CDFM scaling function 
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TABLE I: Values of energies e, angles 9, the coefficient obtained by fitting the magnitude 
of the QE peak, energy shifts e^j^j^ and e^jj^, and transferred momenta q^x^ for the cases of 
inclusive electron scattering cross sections considered. Energies are in MeV, angles are in degrees 
and momenta are in MeV/c. 
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are close to the RFG model results. On the other hand, the results obtained with the use 
of asymmetric CDFM and LFD scaling functions are quite different from those in the RFG 
model, but are close to the predictions of the ERFG model l^^- The basic difference 
from the ERFG model result is observed in the tail extended to small muon energy values, 
where the ERFG model gives more strength. 



V. CONCLUSIONS 



The results of the present work can be summarized as follows: 

i) In Ref. [l^ we extended the CDFM description of the quasielastic ^/''-scaling function 
from 3| by expressing it explicitly and equivalently by means of both density and nucleon 
momentum distributions. In 3, 3| our results on f^^{ip') were obtained on the basis of 
the experimental data on the charge densities for a wide range of nuclei. In the present 
work we extended our approach to consider the scaling function f^^{ip') for medium and 
heavy nuclei with Z ^ N ioT which the proton and neutron densities are not similar. In 
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FIG. 16: The cross section of quasielastic charge-changing reaction [(a), (c)-(f)] and of 

(V^jfi^) reaction (b) on ^^C for e = 1,1.5 and 2 GeV using QE-scahng functions in CDFM (thin 
sohd line: with ci = 0.63; thin dashed line: with ci = 0.72). The results using QE-scahng functions 
in LFD (thick solid line: with ci = 0.63; thick dashed line: with ci = 0.72) are presented in (b) 



and (f). The RFG model result and ERFG result 
lines, respectively. 



i],Q 



are shown by dotted and dash-dotted 
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this case f^^i^'ip') is a sum of the proton and neutron scahng functions calculated by means 
of the proton and neutron densities obtained from nonrelativistic self-consistent mean-field 



calculations. This concerns calculations, as examples, of nuclei 



11 



ike ^^"^Au, ^^Kr, as well as 



12| shows superscaling for 



^^Ni and ^^^Sn [15]. The comparison with the data from 
negative values of the QE if)' including ip' < —1, whereas in the RFG model f{ip') = for 
ip' <-l (see Fig.d). 

ii) We introduce the asymmetry in the CDFM QE scaling function using the fact that the 
maximum value of f{ip') in RFG model is 3/4 while the empirical scaling function reaches 
values smaller than 0.6. In relation with this and the normalization, we parametrize the 
RFG scaling function for V^' > 0, thus simulating the role of all the effects which lead to 
asymmetry and imposing this to the CDFM QE scaling function. In this way, simulating 
phenomenologically the effects which violate the symmetry of f'^^iip') for ip' > including 
the role of the FSI, one can obtain in the CDFM a reasonable agreement of f^^{ip') with 
the empirical data also for positive values of ip' (Fig. |2)). 

iii) We obtain the QE scaling function also on the basis of calculations of nucleon mo- 
mentum distribution n{k) within an approach based on the light-front dynamics method 



45, 46| which improves that used in Ref. \\M- Here we include in the particle-state part 

n n 

of n{k) not only a contribution of the function /s (as in j4^ and [1J|) but also a contribution 
of the function f2- /a and fo are two of the six scalar components of the deuteron total wave 
function in the LED [4^ |4^ and are the main contributions to the tail of Udik). It can be 
seen in Fig. 01 the reasonable agreement of n{k) in LFD with the ?/cvi/- scaling data 0, 0. 
This result made it possible to obtain a good description of the experimental QE scaling 
function (Fig. at least up to ip' ^ —1.2. 

iv) We extend our analysis within the CDFM and LFD to the A-peak region which is the 
main contribution to the inelastic scattering. Here we emphasize that reasonable agreement 
with the experimental data (Fig. was obtained using the empirical value of the coefficient 
in front of the RFG scaling function (0.54 instead of 0.75) in both CDFM and LFD. Also, 
the parameter i?A used in the Fermi-type density for ^^C (necessary to calculate the weight 
function |-Fa(-R)P and thus the scaling function f^{ip'/^)) has a smaller value (1.565 fm) than 
that (i?=2.42 fm) in the QE case, while the value of the diffuseness parameter 6a remains 
the same as h in the QE case. We note that the use of niFnlk) with the same values of (3 
and of kp {(3 = 0.80, kp = 1.20 fm~^) gives a reasonable agreement with results for both 
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QE- and A-region scaling functions (Figs. 0] and EI). 

v) The QE- and A-region scaling functions obtained in the CDFM and in the LED 
approach are applied to description of experimental data on differential cross sections of 
inclusive electron scattering by ^^C at large energies and transferred momenta fEigs. IUin3|). 
The CDEM results are presented for both almost symmetric {c^^ ~ 0.72) and asymmetric 
{c^^ = 0.63) scaling functions. We observe that there are two regions of the value of qf^p 
in different experiments at which the above mentioned (almost symmetric and asymmetric) 
scaling functions work better. The almost symmetric scaling function leads to results in 
agreement with the data in the region of the QE peak in cases when the transferred mo- 
mentum (qf^^) in the position of maximum of the QE peak (ujf^^) is in the scaling region 
ilfci^ > 450 MeV/c ~ 2kp), while the data are overestimated in cases where qf^ < 450 
MeV/c. The results obtained when asymmetric scaling function (c^'^ = 0.63) is used agree 
with the data in cases when < 450 MeV/c and underestimate them when > 450 
MeV/c in the region close to the QE peak, but differences emerge in the transition region. 
In our opinion, the latter case is preferable because additional effects (apart from QE and 
A-resonance) , e.g. of the meson exchange currents could give additional important contri- 
butions to the inclusive electron cross sections for some specific kinematics and minor for 
others. 

vi) The CDEM and LED scaling functions are applied to calculations of QE charge- 
changing neutrino-nuclei reaction cross sections. We present in Eig. the predicted cross 
sections for the reactions {u^, /i~) and (V^, /i+) on the ^^C nucleus for energies of the incident 
particles from 1 to 2 GeV. Our results are compared with those from the REG model and 



from the EREG model |19|,|23l. The results obtained by using the asymmetric CDEM scaling 
function are close to those of EREG and are quite different from the REG results while the 
almost symmetric CDEM scaling function leads to cross sections which are similar to the 
results of the REG model. 
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